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The  formulation  of  the  classical  basic  equations  for  a  thermoelectric  cooler  from  the  Thomson  relations 
to  the  non-linear  differential  equation  with  Onsager’s  reciprocal  relations  was  performed  to  basically 
study  the  Thomson  effect  in  conjunction  with  the  ideal  equation.  The  ideal  equation  is  obtained  when 
the  Thomson  coefficient  is  assumed  to  be  zero.  The  exact  solutions  derived  for  a  commercial  thermo¬ 
electric  cooler  module  provided  the  temperature  distributions  including  the  Thomson  effect.  The  positive 
Thomson  coefficient  led  to  a  slight  improvement  on  the  performance  of  the  thermoelectric  device  while 
the  negative  Thomson  coefficient  led  to  a  slight  declination  of  the  performance.  The  comparison  between 
the  exact  solutions  and  the  ideal  equation  on  the  cooling  power  and  the  coefficient  of  performance  over  a 
wide  range  of  temperature  differences  showed  close  agreement.  In  conclusion,  the  Thomson  effect  is 
small  for  typical  commercial  thermoelectric  coolers  and  the  ideal  equation  effectively  predicts  the 
performance. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  phenomena  are  useful  and  have  drawn  much 
attention  since  the  discovery  of  the  phenomena  in  the  early 
nineteenth  century.  The  barriers  to  applications  were  low  effi¬ 
ciencies  and  the  availability  of  materials.  In  1821,  Thomas  J.  See- 
beck  discovered  that  an  electromotive  force  or  potential  difference 
could  be  produced  by  a  circuit  made  from  two  dissimilar  wires 
when  one  junction  was  heated.  This  is  called  the  Seebeck  effect.  In 
1834,  Jean  Peltier  discovered  the  reverse  process  that  the  passage 
of  an  electric  current  through  a  thermocouple  produces  heating  or 
cooling  depended  on  its  direction  [1].  This  is  called  the  Peltier  ef¬ 
fect  (or  Peltier  cooling ).  In  1854,  William  Thomson  discovered  that 
if  a  temperature  difference  exists  between  any  two  points  of  a 
current-carrying  conductor,  heat  is  either  absorbed  or  liberated 
depending  on  the  direction  of  current  and  material  [2].  This  is 
called  the  Thomson  effect  (or  Thomson  heat).  These  three  effects  are 
called  the  thermoelectric  effects.  Thomson  also  developed  impor¬ 
tant  relationships  between  the  above  three  effects  with  the 
reciprocal  relations  of  the  kinetic  coefficients  (or  sometimes  the 
so-called  symmetry  of  the  kinetic  coefficients)  under  a  peculiar 
assumption  that  the  thermodynamically  reversible  and 
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irreversible  processes  in  thermoelectricity  are  separable  [2].  The 
relationship  developed  is  called  the  Thomson  (or  Kelvin)  relations. 
The  necessity  for  the  assumption  remained  an  objection  to  the 
theory  until  the  advent  of  the  new  thermodynamics  [3-6].  The 
relationship  was  later  completely  confirmed  by  experiments  to  be 
essentially  a  consequence  of  the  Onsager’s  Reciprocal  Principle  in 
1931  [3]. 

As  supported  by  Onsager’s  Principle,  the  Thomson  relations 
provide  a  simple  expression  for  Peltier  cooling,  which  is  the 
product  of  the  Seebeck  coefficient,  the  temperature  at  the  junction, 
and  the  current.  This  Peltier  cooling  is  the  principal  thermoelectric 
cooling  mechanism.  There  are  two  counteracting  phenomena, 
which  are  the  Joule  heating  and  the  thermal  conduction.  The  net 
cooling  power  is  the  Peltier  cooling  minus  these  two  effects. 
Actually  the  Joule  heating  affects  the  thermal  conduction  and 
consequently  the  Peltier  cooling  is  subtracted  only  by  the  thermal 
conduction.  Expressing  the  net  cooling  power  in  terms  of  the  heat 
flux  vector  ~q  gives  [7] 

~q  =  oitJ-  kV  T  (1) 

This  is  the  most  important  formula,  where  a  is  the  Seebeck 
coefficient,  T  the  absolute  temperature,  j  the  current  density  and  k 
the  thermal  conductivity.  Equation  (1)  appears  simple,  but  the 
interpretation  was  somewhat  formidable  because  the  aT  in  the 
equation  was  obtained  by  conjecture  and  the  reciprocal  relations. 
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Nomenclature 

T 

absolute  temperature  (I<) 

Tc 

low  junction  temperature  (I<) 

A 

cross-sectional  area  of  thermoelement  (m2) 

Th 

high  junction  temperature  (I<) 

COP 

the  coefficient  of  performance,  dimensionless 

h 

low  junction  temperature  (I<) 

cp 

specific  heat  capacity  (J/kg  I<) 

t2 

high  junction  temperature  (I<) 

1 

electric  field  vector  (V/m) 

AT 

temperature  difference  (I<),  T2  -  T\ 

I 

electric  current  (A) 

V 

volume  (m3) 

W 

work  per  unit  time  W  =  alAT  + 12 R 

j 

electric  current  density  vector  (A/m2) 

X 

distance  of  thermoelement  leg  (m) 

j 

electric  current  density  (A/m2) 

arbitrary  quantities 

Ji 

I< 

arbitrary  quantities 
thermal  conductance  (W/K) 

Xi 

arbitrary  quantities:  dXi/dt 

L 

length  of  thermoelement  (m) 

C,vppP  svmhnfc 

k 

thermal  conductivity  (W/m  I<) 

a 

Seebeck  coefficient  (V/K) 

heat  flux  vector  (W/m2) 

P 

dimensionless  number  =  fT2(da/dT)/(Ak/L) 

Q 

heat  generation  (W/m3) 

T 

dimensionless  number  7  =  72ft/(A/<AT/L) 

Qc 

cooling  power  (W) 

0 

the  ratio  of  AT  to  T2 

Qi 

cooling  power  (W) 

P 

density  (kg/m3) 

R 

electrical  resistivity  (Q) 

P 

electric  resistivity  (Q  m) 

Rij 

kinetic  coefficients 

T 

Thomson  coefficient  (V/I<) 

%en 

entropy  generation  per  unit  volume  (J/m3  I<) 

6 

dimensionless  temperature  0  =  (T  -  Ti)/(T2  -  Ti) 

Sgen 

entropy  generation  per  unit  time  and  unit  volume 

? 

dimensionless  distance  ?  =  x/L 

CJ/m3  K) 

@1 

dimensionless  cooling  power  ©!  =  Qi/(A/<AT/L) 

Sgen 

entropy  generation  (J/I<) 

n 

dimensionless  Peltier  cooling  77  =  a(T1)//(A/c/L) 

Sgen 

entropy  generation  per  unit  time  (W/I<) 

w 

dimensionless  work  W  =  W/(AkAT/L ) 

In  fact,  the  three  thermoelectric  effects  that  are  known  to  be 
reversible  are  effectively  combined  into  the  aT  [8-10].  The  second 
term  gives  the  thermal  conduction  that  is  obviously  irreversible. 
A  steady-state  heat  diffusion  equation  is  given  by  [7] 

V  (kv  rj+jVr^.v  t  =  o  (2) 

where  p  is  the  electrical  resistivity.  Equation  (2)  is  a  non-linear 
differential  equation.  The  Thomson  coefficient  is  given  by 


In  order  to  obtain  the  net  cooling  power,  Equation  (2)  is  first 
solved  for  the  thermal  conduction  including  the  Joule  heating  and 
then  Equation  (1)  can  give  the  cooling  power  at  the  junction.  It  is 
often  assumed  [8-10,16,17]  that  the  Thomson  effect  was  absent  or 
negligible  or,  in  other  words,  a  was  independent  of  temperature.  In 
such  a  case,  Equation  (2)  becomes  an  ordinary  differential  equation 
and  can  be  easily  solved  for  the  thermal  conduction  and  then 
Equation  (1)  yields  a  well-known  equation  as 

Qc  =  aTcI-^I2R-K(Th-Tc)  (4) 

where  Qc  is  the  cooling  power,  I  the  current,  Tc  the  low  junction 
temperature,  Th  the  high  junction  temperature,  R  the  electrical 
resistance,  and  I<  the  thermal  conductance.  Equation  (4)  is  called 
the  ideal  equation  and  has  been  widely  used  in  science  and  in¬ 
dustry  [8,10,15,16]  often  with  a  fair  agreement  with  the  experi¬ 
ments  [9,13]. 

The  third  term  in  Equation  (2)  was  regarded  as  the  Thomson 
effect.  The  influence  of  the  Thomson  effect  was  discussed  in  the 
literature  [11-14]  by  numerically  solving  Equation  (2),  in  which 
the  Seebeck  coefficient  a  was  dependent  of  temperature.  The 


exact  solutions  of  Equation  (2)  assuming  that  the  Thomson  co¬ 
efficient  t  is  constant  indicated  that  the  Thomson  effect  would 
improve  the  performance  significantly  [12,13].  More  realistic 
work  on  the  Thomson  coefficient  assuming  that  da/dT  is  a  linear 
function  of  temperature  was  numerically  conducted  using  a  finite 
difference  method  and  also  experimentally  compared  with  a 
commercial  module  showing  slight  improvement  on  the  perfor¬ 
mance  [13].  Numerical  simulations  on  the  Thomson  effect  using 
CFD  software  ANSYS  for  miniature  thermoelectric  coolers 
revealed  that  the  cooling  power  can  be  improved  by  a  factor  of 
5—7%  [14]  by  including  the  Thomson  effect.  There  is  an  experi¬ 
mental  work  [18]  showing  that  the  Thomson  effect  is  responsible 
for  reduction  of  the  figure  of  merit  on  a  thermoelectric  generator 
with  an  increase  of  temperature  difference,  whereas  the  re¬ 
searchers  anticipated  enhancement  of  the  figure  of  merit  on  a 
thermoelectric  cooler.  There  were  attempts  [19,20]  to  improve 
the  performance  of  thermoelectric  devices  by  optimizing  the 
temperature  dependence  of  the  material  properties  wherein  the 
Thomson  effect  was  taken  into  account  although  the  improve¬ 
ment  could  not  be  demonstrated  by  the  measurements.  There  is  a 
rigorous  numerical  analysis  [21]  for  emerging  materials  and  large 
temperature  gradients  in  thermoelectric  generators  how  well  the 
ideal  equation  predicts  the  exact  solutions  including  the  Thom¬ 
son  effect.  It  is  analytically  proven  [21]  that  the  ideal  equation 
produces  the  exact  solution  output  power  and  efficiency  despite 
its  limiting  assumptions  if  an  integral-averaged  Seebeck  coeffi¬ 
cient  is  used. 

From  the  review  of  the  above  experimental  and  theoretical 
work,  the  present  author  feels  that  further  realistic  studies  are 
required  to  clarify  the  Thomson  effect  related  to  the  ideal  equation 
in  temperature  distributions,  particularly  in  the  positive  and 
negative  Thomson  coefficients.  The  present  work  studies  the 
Thomson  effect  in  conjunction  with  the  ideal  equation  from  the 
formulation  of  the  basic  equations  to  the  numerical  simulations  in 
light  of  a  real  commercial  thermoelectric  module. 
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2.  Formulation  of  basic  equations 


2.2.  Basic  equations 


2  A.  Onsager’s  reciprocal  relations 

The  second  law  of  thermodynamics  with  no  mass  transfer 
(AS  =  0)  in  an  isotropic  substance  provides  an  expression  for  the 
entropy  generation. 

2 

Sgen  =  —  J  ~y~  >  0  (5) 

1 

Suppose  that  the  derivatives  of  the  entropy  generation  sgen  per 
unit  volume  with  respect  to  arbitrary  quantities  x2  are  quantities  Ji 
[4,7]  as 


Let  us  consider  a  non-uniformly  heated  thermoelectric  material. 
For  an  isotropic  substance,  the  continuity  equation  for  a  constant 
current  gives 

V  T  =  0  (16) 

The  electric  field  I?  js  affected  by  the  current  density  7  and  the 
temperature  gradient  V  T.  The  coefficients  are  known  according  to 
the  Ohm’s  law  and  the  Seebeck  effect  [7].  The  field  is  then 
expressed  as 

If  =  Jp  +  aV  T  (17) 


d5gen  _  r 
dXi  Ji 


where  p  is  the  electrical  resistivity  and  a  is  the  Seebeck  coefficient. 
Solving  for  the  current  density  gives 


At  maximum  sgen,  Ji  is  zero.  Accordingly,  at  the  state  close  to 
equilibrium  (or  maximum),  J*  is  small  [4].  The  entropy  generation 
per  unit  time  and  unit  volume  is 


•  _  dSgen  _  d*j  j_ 

gen  ~  dt  ~  dtJi 

We  want  to  express  dxjdt  as  X*.  Then,  we  have 


(7) 


-•Sgen  —  XJi  (8) 

The  X/  are  usually  a  function  of  Ji.  Onsager  [3]  stated  that  if  Ji 
were  completely  independent,  we  should  have  the  relation, 
expanding X/  in  the  powers  of  Ji  and  taking  only  linear  terms  [7].  The 
smallness  of  Ji  in  practice  allows  the  linear  relations  being  sufficient. 


%  = 

j 


(9) 


7  =  ’(£-«VT)  (18) 

The  heat  flux  qMs  also  affected  by  the  both  field  E  and  tem¬ 
perature  gradient  V  T.  However,  the  coefficients  were  not  readily 
attainable  at  that  time.  Thomson  in  1854  arrived  at  the  relationship 
assuming  that  thermoelectric  phenomena  and  thermal  conduction 
are  independent.  Later,  Onsager  [3]  supported  the  relationship  by 
presenting  the  reciprocal  principle  which  was  experimentally 
proved  but  failed  the  phenomenological  proof  [6].  Here  we  derive 
the  Thomson  relationship  using  Onsager’s  principle.  The  heat  flux 
with  arbitrary  coefficients  is 


q  =C,E  -C2VT  (19) 

The  general  heat  diffusion  equation  is 


where  Ry  are  called  the  kinetic  coefficients.  Hence,  we  have 

-Sgen  =  CO) 

i 

It  is  necessary  to  choose  the  quantities  X2-  in  some  manner,  and 
then  to  determine  the  Ji.  The  X;  and  Ji  are  conveniently  determined 
simply  by  means  of  the  formula  for  the  rate  of  change  of  the  total 


entropy  generation. 

Sgen  =  -  f  Y)XJtdV  >  0  (11) 

For  two  terms,  we  have 

Sgen  =  -  J[X,h  +X2j2]dV  >  0  (12) 

And 

Xi  =  R11J1+R12I2  (13) 

%2  =  *2l/l  +R22J2  (14) 

From  these  equations,  one  can  assert  that  the  kinetic  coefficients 
are  symmetrical  with  respect  to  the  suffixes  1  and  2. 

Ru  =  *21  (15) 

which  is  called  the  reciprocal  relations. 


-V-q+q  =  pcp—  (20) 

For  steady  state,  we  have 

-V  ~q  +q  =  0  (21) 

where  q  is  expressed  by  [7] 

q  =  if- 7  =  J2p  +  7*aV  T  (22) 

Expressing  Equation  (5)  with  Equation  (21)  in  a  manner  that  the 
heat  flux  gradient  V  •  ~q  has  a  minus  sign,  the  entropy  generation 
per  unit  time  will  be 

2  . 

Sgen  =  -jd~Y  =  Jj(-v  A  +  q)dv  (23) 

1 

Using  Equation  (22),  we  have 

Sgen  =  J  -  J  l^dV  (24) 

The  second  term  is  integrated  by  parts,  using  the  divergence 
theorem  and  noting  that  the  fully  transported  heat  does  not  pro¬ 
duce  entropy  since  the  volume  term  cancels  out  the  surface 
contribution  [7].  Then,  we  have 
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Note  that  the  Joule  heating  and  the  thermal  conduction  in  the 
equation  are  indeed  irreversible.  This  equation  satisfies  the 
requirement  for  Onsager’s  reciprocal  relations  as  shown  in  Equation 
(12).  Substituting  Equations  (22)  and  (30)  in  Equation  (21)  yields 


V  (kv  T)  +J'2p -r“|7'V  T  =  0  (32) 

The  Thomson  coefficient  t,  originally  obtained  from  the  Thom¬ 
son  relations,  is  written 


(33) 


In  Equation  (32),  the  first  term  is  the  thermal  conduction,  the 
second  term  is  the  Joule  heating,  and  the  third  term  is  the  Thomson 
heat.  Note  that  if  the  Seebeck  coefficient  a  is  independent  of  tem¬ 
perature,  the  Thomson  coefficient  t  is  zero  and  then  the  Thomson 
heat  is  absent. 


If  we  take  j  and  ~q  as  X\  and  X2,  then  the  corresponding 
quantities  J\  and  J2  are  the  components  of  the  vectors  -E/E  and 
V  T /T2  [7].  Accordingly,  as  shown  in  Equations  (13)  and  (14)  with 
Equation  (19),  we  have 


T  =  -c2r2^T  (27) 

where  the  reciprocal  relations  in  Equation  (15)  gives 


2.3.  Exact  solutions 


Let  us  consider  one  of  p-  and  n-type  thermoelements  as  shown 
in  Fig.  1,  knowing  that  p-  and  n-type  thermoelements  produce  the 
same  results  if  the  materials  and  dimensions  are  assumed  to  be 
similar.  Equation  (32)  for  one  dimensional  analysis  at  steady  state 
gives  the  non-linear  differential  governing  equation  as 


dx2  Ak  dx  A2/< 


(34) 


where  A  is  the  cross-sectional  area  of  the  thermoelement.  We  want 
to  make  it  dimensionless.  The  boundary  conditions  will  be  T(0)  =  T\ 
and  T(L)  =  T2.  Then,  Let 

d  =  T^andZ  =  I  W 


QT  =  -aT2  (28) 

P 

Thus,  we  have  Cj  =  aT/p.  Inserting  this  and  Equation  (17)  into  (27) 
gives 


where  L  is  the  element  length.  Then,  the  dimensionless  differential 
equation  is 


d  2_d 
d'2 


@((6  —  ~i)4>  +  \)  gi'+y 


o 


(36) 


t  =  «t7-(c2-^V7  (29) 

We  introduce  the  Wiedemann-Franz  law  that  is  L0  =  kp/T , 
where  L0  is  the  Lorentz  number  (constant).  Expressing  k  as  L0T/p,  we 
find  that  the  coefficient  (C2  -a2T/p)  is  nothing  more  than  the 
thermal  conductivity  k.  Finally,  the  heat  flow  density  vector  (heat 
flux)  is  expressed  as 

~q  =  oltJ  -  kV  T  (30) 

This  is  the  most  essential  equation  for  thermoelectric  phe¬ 
nomena.  The  second  term  pertains  to  the  thermal  conduction  and 
the  term  of  interest  is  the  first  term,  which  gives  the  thermoelectric 
effects:  directly  the  Peltier  cooling  but  indirectly  the  Seebeck  effect 
and  the  Thomson  heat  [8-10].  Inserting  Equations  (17)  and  (30) 
into  Equation  (25)  gives 

Sgen  =  j(^f  +  ^rP~yV>0  (31) 


where  the  dimensionless  numbers  are  defined  as  follows.  The 
boundary  condition  will  be  that  6(0)  =  0  and  6(  1)  =  1.  Equation 
(36)  was  formulated  on  the  basis  of  the  high  junction  temperature 
T2  to  meet  the  commercial  products. 


P  = 


ITJtL 

Ak  — 


(37) 


where  is  the  ratio  of  the  Thomson  heat  to  the  thermal  conduction. 
Note  that  is  not  a  function  of  AT. 


T  = 


I2R 
Ak- — 


(38) 


where  y  is  the  ratio  of  the  Joule  heating  to  the  thermal  conduction. 
AT 

0  =  (39) 


The  entropy  generation  per  unit  time  for  the  irreversible 
processes  is  indeed  greater  than  zero  since  p  and  k  are  positive. 


where  0  is  the  ratio  of  temperature  difference  to  the  high  junction 
temperature.  The  temperature  difference  is 
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AT  =  T2-  Ta  (40) 

where  AT  is  the  high  junction  temperature  T2  minus  the  low 
junction  temperature  T\.  Therefore,  T\  is  a  function  of  AT  since  T2  is 
fixed.  The  cooling  power  at  the  cold  junction  using  Equation  (30)  is 


at  the  average  temperature  of  AT  [15,16]  as  a  result  of  the  internal 
phenomena  being  imposed  on  the  surface  phenomena. 


3.  Results  and  discussion 


Qi 


«(Ti)Ti /  + 


,  .dT 

-  kA-j- 
dx 


x=0 


(41) 


where  the  first  term  is  the  Peltier  cooling  and  the  second  term  is  the 
thermal  conduction.  It  has  been  customary  in  the  literature  for  the 
exact  solution  wherein  the  Seebeck  coefficient  a  is  evaluated  at  the 
cold  junction  temperature  T\.  The  dimensionless  cooling  power  is 
expressed 


©i 


(42) 


where  6^  is  the  dimensionless  cooling  power,  which  is 


1  AkAT/L 


(43) 


and  the  dimensionless  Peltier  cooling  is 


tt  =  a(T-,)I 
Ak/L 


(44) 


The  work  per  unit  time  gives  W  =  alAT  +  /2ft  [10].  Then  the 
dimensionless  work  per  unit  time  W  is  expressed  by 


W  =  77  +  7 


(45) 


where  W  =  W/ (AkAT/L).  Then,  the  coefficient  of  performance 
(COP)  for  the  thermoelectric  cooler  is  determined  as  the  cooling 
power  over  the  work  as 

COP  =  |l  (46) 

Equation  (36)  can  be  exactly  solved  for  the  temperature  distri¬ 
butions  conveniently  with  mathematical  software  Mathcad  and 
then  the  cooling  power  of  Equation  (41 )  can  be  obtained,  which  are 
the  exact  solutions  including  the  Thomson  effect. 


2.4.  The  ideal  equation 

According  to  the  assumption  made  by  the  both  Thomson  re¬ 
lations  and  Onsager’s  reciprocal  relations,  the  thermoelectric  ef¬ 
fects  and  the  thermal  conduction  in  Equation  (30)  are  completely 
independent,  which  implies  that  each  term  may  be  separately  dealt 
with.  In  fact,  separately  dealing  with  the  reversible  processes  of  the 
three  thermoelectric  effects  yielded  the  Peltier  cooling  of  the  aTI  in 
the  equation  [8-10].  However,  under  the  assumption  that  the 
Thomson  coefficient  is  negligible  or  the  Seebeck  coefficient  is  in¬ 
dependent  of  temperature,  Equation  (34)  easily  provides  the  exact 
solution  of  the  temperature  distribution.  Then,  Equation  (41 )  with 
the  temperature  distribution  for  the  cooling  power  at  the  cold 
junction  reduces  to 

1  Ak 

Qi  =  a(Tavg)T1I-^I2R--j-(T2-T,)  (47) 

which  appears  simple  but  is  a  robust  equation  with  a  usually  good 
agreement  with  experiments  and  with  comprehensive  applications 
in  science  and  industry.  This  is  here  called  the  ideal  equation.  The 
ideal  equation  assumes  that  the  Seebeck  coefficient  a  be  evaluated 


In  order  to  realistically  study  the  Thomson  effect  on  the  tem¬ 
perature  distributions  with  the  temperature  dependent  Seebeck 
coefficient,  a  real  commercial  module  of  Laird  CP10-127-05  was 
chosen.  Both  the  temperature  dependent  properties  and  the  di¬ 
mensions  for  the  module  were  provided  by  the  manufacturer, 
which  are  shown  in  Table  1  and  the  temperature  dependent  See¬ 
beck  coefficient  is  graphed  in  Fig.  2. 

To  reveal  the  Thomson  effect,  only  the  Seebeck  coefficients  are 
considered  to  be  dependent  of  the  temperature  while  the  electrical 
resistivity  and  the  thermal  conductivity  are  constant  as  shown  in 
Table  1.  In  Fig.  2,  the  Seebeck  coefficient  increases  with  increasing 
the  temperature  up  to  about  350  K  and  then  decreases  with 
increasing  the  temperature,  so  that  there  is  a  peak  value.  For  the 
feasible  solution  with  the  non-linear  differential  equation,  the 
linear  curve  fits  for  the  Seebeck  coefficient  were  used  in  the  present 
work.  In  Fig.  2,  the  linear  curve  fit  of  the  first  part  for  a  range  from 
250  K  to  350  K  was  used  in  the  present  work,  while  the  other  linear 
curve  fit  of  the  second  part  for  the  range  from  350  K  to  450  K  was 
used  separately.  The  linear  curve  fits  of  the  two  parts  are  shown  in 
Table  1.  Some  reports  [11-14]  illustrated  used  only  the  mono- 
tonically  increasing  Seebeck  coefficients.  Therefore,  the  negative 
slope  of  the  Seebeck  coefficient  was  not  presented.  As  mentioned 
previously,  the  absolute  values  of  the  temperature  dependent 
Seebeck  coefficients  for  p-  and  n-type  elements  are  assumed  to  be 
the  same,  but  the  sign  of  n-type  element’s  coefficient  is  negative 
while  the  sign  of  p-type  element’s  coefficient  is  positive.  In  reality, 
the  Seebeck  coefficients  in  the  p-  and  n-type  elements  are  different 
as  the  absolute  values  vary  with  temperature.  Fig.  2  would  be  the 
arithmetic  average  of  the  two.  By  doping  technique  (intentionally 
introducing  impurities  into  a  pure  semiconductor),  n-type  semi¬ 
conductor  carries  predominantly  free  electrons  moving  to  the 
opposite  direction  of  current  applied  (which  causes  the  negative 
sign  of  the  Seebeck  coefficient),  while  p-type  semiconductor  carries 
free  holes  in  the  same  manner  but  moving  to  the  same  direction  of 
current  applied  (which  causes  the  positive  sign).  Likewise,  p-  and 
n-type  thermoelements  can  be  arranged  electrically  in  series  and 
thermally  in  parallel  to  form  a  module.  The  number  of  electrons  and 
holes  can  be  controlled  through  the  introduction  of  impurities 
(doping),  which  are  called  donors  or  acceptors. 

The  dimensionless  differential  equation,  Equation  (36),  was 
developed  on  the  basis  of  a  fixed  high  junction  temperature  T2  so 
that  the  low  junction  temperature  T\  decreases  as  AT  increases.  0  in 
Equation  (37)  is  the  dimensionless  number  indicating  approxi¬ 
mately  the  ratio  of  the  Thomson  heat  to  the  thermal  conduction.  0 
is  a  function  only  of  the  slope  da/dT and  the  current  7  Equation  (36) 


Table  1 

Maximum  values,  dimensions,  and  the  properties  for  a  commercial  module. 


TEC  Module  (Laird  CPI 0-1 27-05) 

At  25  °C 

ATmax  (°C) 

67 

Tnax  (A) 

3.9 

Qmax  (W) 

34.3 

L  (mm),  element  length 

1.253 

A  (mm2),  cross-sectional  area 

1 

#  of  thermocouples 

127 

a  (T),  (pV/K) 

0.2068T  +  138.78,  250  K  <  T  <  350  K 

a  CO,  (pV/K) 

— 0.2144T  +  281.02,  350  K  <  T  <  450  K 

p  (Q  cm)  at  27  °C 

1.01  x  10"3 

k  (W/cmK)  at  27  °C 

1.51  x  10"2 

Module  size 

30  x  30  x  3.2  mm 

66 
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with  the  given  boundary  conditions  was  solved  for  the  temperature 
distributions  conveniently  using  mathematical  software  Mathcad. 

A  typical  value  of  ^  =  0.2  was  obtained  for  the  commercial 
product  at  T2  =  298  K  and  /max  =  3.9  A.  And  a  hypothetical  value  of 
1.0  that  is  fivefold  larger  than  the  typical  value  of  0.2  was  used  in 
order  to  closely  examine  the  Thomson  effect,  y  in  Equation  (38)  is 
the  dimensionless  number  indicating  approximately  the  ratio  of 
the  Joule  heating  to  the  thermal  conduction. 

Fig.  3(a)  and  (b)  show  the  temperature  distributions  at 
T2  =  298  K  and  /  =  3.9  A  as  a  function  of  (3  for  AT  =  10  K  (y  =  15.97) 


l  =  x/L 


Fig.  3.  (a)  AT  =  10  K  (y  =  15.97),  and  (b)  AT  =  50  K  (y  =  3.19),  dimensionless  tem¬ 
perature  6  vs.  dimensionless  distance  £  as  a  function  of  0.  The  operating  conditions  are 
T2  -  298  K  and  /  =  3.9  A. 


and  AT  =  50  K  (y  =  3.19).  In  Fig.  3(a),  the  Joule  heating  appears 
dominant  with  a  high  value  of  y  =  15.97,  while  in  Fig.  3(b)  the 
thermal  conduction  is  approximately  balanced  with  the  Joule 
heating  with  a  low  value  of  y  =  3.19.  For  the  commercial  cases  in 
the  both  figures  with  (3  =  0.2,  the  Thomson  effects  appear  small 
compared  to  the  ideal  equation  at  0  =  0.  The  hypothetical  cases  at 
=  1.0  obviously  provide  salient  examples  for  the  Thomson  effect. 
Since  the  temperature  gradient  at  ?  =  0  counteracts  the  Peltier 
cooling,  the  slightly  lower  temperature  distribution  near  zero  acts 
as  improving  the  net  cooling  power.  The  Thomson  heat  is  effec¬ 
tively  described  by  the  product  of  0  and  d0/d?.  The  Thomson  heat 
acts  as  absorbing  heat  when  the  fidO/dt,  is  positive,  while  it  acts  as 
liberating  heat  when  the  /?d0/d?  is  negative,  which  are  shown  in 
Fig.  3(a).  In  the  first  half  of  the  temperature  distribution,  both  and 
d$/d?  are  positive,  so  that  the  product  is  positive.  It  is  known  that 
the  moving  charged  electrons  or  holes  (see  Fig.  1 )  transport  not  only 
the  electric  energy  but  also  the  thermal  energy  as  absorbing  and 
liberating  depending  on  the  sign  of  (3dd/d%  along  the  thermoele¬ 
ment.  The  temperature  distribution  is  slightly  shifted  to  the  right 
giving  the  lower  temperature  distribution  near  f  =  0,  which  results 
in  the  improved  cooling  power  at  the  cold  junction.  It  is  interesting 
to  note  that  the  product  of  the  current  I  and  the  Seebeck  coefficient 
a  in  Equation  (37)  is  always  positive,  therefore  making  no  contri¬ 
bution  to  the  final  sign  of  the  Thomson  heat.  If  the  current  is  pos¬ 
itive  in  a  p-type  element,  the  sign  of  a  is  positive.  If  the  current  is 
negative  in  an  n-type  element,  the  sign  of  a  is  negative  (see  Fig.  1). 
Therefore,  the  product  is  always  positive. 

Now  consider  the  linear  curve  fit  of  the  second  part  (for  the 
range  from  350  K  to  450  K)  of  the  temperature-dependent  Seebeck 
coefficient  in  Fig.  2,  where  the  (3  in  Equation  (37)  is  negative. 
Fig.  4(a)  and  (b)  show  the  temperature  distributions  at  T2  =  420  I< 
and  /  =  3.9  A  as  a  function  of  /3  for  AT  =  10  K  (y  =  15.97)  and 
AT  =  50  K  (y  =  3.19).  For  the  commercial  cases  at  (3  =  -0.29  in  the 
both  figures,  the  Thomson  effect  again  appears  small  compared  to 
the  ideal  equation  at  0  =  0.  The  hypothetical  cases  at  /3  =  -1.46 
obviously  provide  salient  examples  for  the  Thomson  effect.  Note 
that  the  temperature  distribution  is  slightly  shifted  to  the  left  giving 
the  higher  temperature  gradient  near  f  =  0,  which  results  in  the 
reduced  cooling  power  at  the  cold  junction.  This  is  opposite  of  the 
cases  at  T2  =  298  I<  shown  in  Fig.  3(a)  and  (b).  In  Fig.  4(a),  it  is  seen 
that  the  first  and  second  half  segments  of  the  thermoelement 
exhibit  the  heat  liberation  and  absorption,  while  the  moving 
charged  carriers  (electrons  or  holes)  are  moving  from  left  to  right. 
Actually  the  moving  carriers  are  resultantly  at  a  speed  from  left  to 
right,  while  numerous  carriers  are  still  moving  in  random  di¬ 
rections,  which  can  transport  the  thermal  energy  from  right  to  left 
without  problems. 

Equation  (41 )  represents  the  net  cooling  power  which  is  calcu¬ 
lated  by  subtracting  the  thermal  conduction  from  the  Peltier 
cooling.  Each  term  of  the  equation  is  separately  presented  in 
Fig.  5(a)  and  (b)  in  order  to  examine  the  fraction  of  the  thermal 
conduction  to  the  Peltier  cooling  as  well  as  the  Thomson  effect  in 
terms  of  the  heat  transfer  rate.  In  Fig.  5(a)  for  AT  =  10  K,  the  per¬ 
centage  of  the  thermal  conduction  to  the  Peltier  cooling  at  I  =  3.9  A 
is  46.2%,  while,  in  Fig.  5(b)  for  AT  =  50  K,  the  percentage  is  81.6%. 
This  indicates  that  the  significant  portion  of  heat  generated  by  the 
Peltier  cooling  is  wasted  surprisingly  by  the  thermal  conduction 
and  this  waste  robustly  increases  as  AT  increases.  Note  that  the  net 
cooling  power  is  zero  at  I  =  1.5  A.  No  gain  in  the  cooling  power  is 
seen  until  after  I  =  1.5  A.  It  would  be  desirable  from  the  above 
cognition  to  develop  a  thermoelectric  material  that  has  a  low 
thermal  conductivity  to  minimize  the  waste  heat  by  the  thermal 
conduction.  However,  lowering  the  thermal  conductivity  without 
affecting  the  electrical  conductivity  is  challenging  according  to  the 
Wiedemann— Franz  law. 
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£,  =  x/L 


£,  =  x/L 

Fig.  4.  (a)  AT  =  10  K  (y  =  15.97),  and  (b)  AT  =  50  K  (y  =  3.19),  dimensionless  tem¬ 
perature  6  vs.  dimensionless  distance  £  as  a  function  of  f3.  The  operating  conditions  are 
T2  =  420  K  and  I  =  3.9  A. 


The  net  cooling  power  Qi  in  Equation  (41 )  and  the  coefficient  of 
performance  (COP)  in  Equation  (46)  along  the  current  at  both 
T2  =  298  K  (thin  lines)  and  420  K  (thick  lines)  for  AT  =  10  I<  are 
presented  in  Fig.  6.  The  maximum  current  for  this  commercial 
module  is  3.9  A  as  shown  in  Table  1,  so  the  operating  current  barely 
exceeds  the  maximum  current  in  practice.  The  slight  improvement 
by  the  Thomson  effect  on  Qi  at  298  K  is  also  reflected  on  the  cor¬ 
responding  COP  wherein  the  effect  appears  even  smaller.  Of  inter¬ 
est  is  that  Qi  at  T2  =  420  K  makes  in  the  opposite  way.  The  Thomson 
effect  improves  Q.i  at  298  K  while  it  deteriorates  at  420  K. 

As  A T  increases,  the  both  cooling  power  and  COP  decrease, 
which  is  shown  in  Fig.  7.  There  is  no  appreciable  change  on  the 
Thomson  effect  with  increasing  AT. 

The  net  cooling  powers  at  T2  =  298  K  and  420  I<  as  a  function  of 
AT  are  presented  in  Fig.  8.  At  a  glance  the  exact  solutions  (jS  ¥=  0) 
and  the  ideal  equation  (/?  =  0)  are  very  close  to  each  other. 
Furthermore,  the  two  at  298  K  are  in  good  agreement  with  the 
available  data  of  commercial  products  which  are  usually  empirical 
data.  It  should  be  pointed  out  that  the  gaps  between  the  Thomson 
effect  and  the  ideal  equation  at  both  T2  =  298  K  and  420  K  remain 
unchanged,  which  is  interpreted  to  be  good  for  the  predictability  by 
the  ideal  equation.  It  is  postulated  that  the  ideal  equation  is  very 
effective  for  prediction  of  the  performance  of  thermoelectric 
coolers,  if  the  averaged  Seebeck  coefficient  is  used. 

In  order  to  elucidate  the  agreement  between  the  exact  solutions, 
the  ideal  equation,  and  the  numerical  simulations,  the  numerical 
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Fig.  5.  (a)  AT  =  10  K  (y  =  15.97),  and  (b)  AT  =  50  K  (y  =  3.19),  heat  transfer  rate  vs. 
current  at  T2  =  298  K. 
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Fig.  6.  Cooling  power  and  COP  vs.  current  at  T2  =  298  I<  and  420  K  for  AT  =  10  K 
0  =  -0.29  at  T2  -  420  K  and  0  =  0.2  at  T2  =  298  K. 
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Fig.  7.  Cooling  power  and  COP  vs.  current  at  T2  =  298  K  and  T2  =  420  K  for  AT  =  50  K 
0  -  -0.29  at  T2  =  420  K  and  §  =  0.2  at  T2  =  298  K. 
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software  ANSYS  was  deployed  with  the  exactly  same  material 
properties,  dimensions,  and  boundary  conditions  as  those  imposed 
on  the  exact  solutions.  Recall  that  the  exact  solutions  were 
analytically  obtained  using  mathematical  software  Mathcad. 

The  exact  solutions  with  Equations  (36)  and  (42)  and  the 
methodology  used  are  justified  by  comparing  them  with  the  ANSYS 
simulations  as  shown  in  Fig.  9.  The  ANSYS  solves  numerically  the 
governing  equations  of  Equations  (17),  (20),  (22),  and  (30)  by  means 
of  the  finite  element  method  with  the  given  boundary  conditions 
[17].  It  is  found  in  Fig.  9  that  the  exact  solutions  are  in  precise 
agreement  with  the  ANSYS  simulations.  Accordingly,  the  gaps  be¬ 
tween  the  exact  solutions  and  the  ideal  equation  implicate  the 
quantities  associated  with  the  assumption  made  in  the  ideal 
equation.  It  is  also  deduced  that  the  assumption  of  the  negligible 
Thomson  coefficient  is  indeed  appropriate  for  either  the  moderate 
currents  or  the  moderate  temperature  differences.  However,  the 
errors  between  the  ideal  equation  and  the  exact  solutions  could 
reach  up  to  8%  for  the  maximum  current  (here  J  =  3.9  A)  and  the 
high  temperature  difference  (AT  =  50  I<).  All  the  procedures  and 
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Fig.  8.  Cooling  power  as  a  function  of  AT  at  T2  =  298  K  and  420  K  for  I  =  3.9  A. 
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Fig.  9.  Comparison  of  the  exact  solution  with  ANSYS  simulations  at  T2  =  298  K. 

conditions  implemented  in  the  exact  solutions  and  the  ideal 
equation  turned  out  to  be  correct. 

4.  Conclusions 

A  positive  Thomson  coefficient  improves  the  performance  of  a 
thermoelectric  cooler  while  a  negative  Thomson  coefficient  re¬ 
duces  the  performance. 

A  substantial  amount  of  heat  absorbed  by  the  Peltier  cooling  is 
wasted  by  the  thermal  conduction,  substantially  increasing  with 
increasing  the  temperature  difference. 

The  exact  solutions  of  the  analytical  differential  equations 
including  the  Thomson  effect  are  in  excellent  agreement  with  the 
ANSYS  numerical  simulations. 

The  ideal  equation  effectively  predicts  the  performances  of 
thermoelectric  coolers  if  the  Seebeck  coefficient  is  correctly  eval¬ 
uated  at  the  average  temperature  of  the  low  and  high  junction 
temperatures. 
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